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ABSTRACT 

Galaxies are believed to be in one-to-one correspondence with simulated dark mat- 
ter subhalos. We use high-resolution N-body simulations of cosmological volumes to 
calculate the statistical properties of subhalo (galaxy) major mergers at high redshift 
(z — 0.6 — 5). We measure the evolution of the galaxy merger rate, finding that it is 
much shallower than the merger rate of dark matter host halos at z > 2.5, but roughly 
parallels that of halos at z < 1.6. We also track the detailed merger histories of indi- 
vidual galaxies and measure the likelihood of multiple mergers per halo or subhalo. We 
examine satellite merger statistics in detail: 15% — 35% of all recently merged galaxies 
are satellites and satellites are twice as likely as centrals to have had a recent major 
merger. Finally, we show how the differing evolution of the merger rates of halos and 
galaxies leads to the evolution of the average satellite occupation per halo, noting that 
for a fixed halo mass, the satellite halo occupation peaks at z ~ 2.5. 

Key Vifords: cosmology:theory - methods:N-body simulations - galaxies: halos - 
galaxies: interactions 



1 INTRODUCTION 



Mergers are key in the hierarchical growth of structure, 
and major galaxy mergers (referred to as mergers hence- 
forth) are thought to play a crucial role in galaxy evolu- 
tion. Specifically, the y are expected to t rigger quasar activity 
(e.g.. ICarlber3ll99d) . starbursts (e.g., iBarnes fc HernauistI 
199l|), and morphological changes (e.g.. 



iGhigna et al.l [ioOCj lFl Simulations are now becoming suffi- 
ciently high in resolution and large in volume to provide 
statistically significant samples (e.g. [ Pe Lucia et al. |2004| : 
iDiemand etloT l2004l : iGao et al.l l2004l : iReed et al.l boosf ). 



199li:7Noguchil 



Toomre fc Toomrd \l97^ . and they are thought to be 
related to Lyman Break Galaxies (LBGs), submillime- 
ter galaxies (SMGs), an d ultra-luminous infrared galaxies 
(ULIRGs) (see reviews bylGiavaliscdl2002l : iBlain et aLll2002l : 
ISanders fc Mirabellll99B . respectively). Observational sam- 
ples of such objects at z > 1 are now becoming large 



This high mass and force resolution is necessary to track 
bound subhalos throughout their orbit in the host halo and 
avoid artificial numerical disruption. This is particularly im- 
portant for tracking the orbits of galaxies, which are ex- 
pected to reside in the dense inner core of subhalos and to 
be more stable to mass stripping than dark matter because 
of dissipative gas dynamics. The correspondence of galax- 
ies with subhalos has been successful in reproducing galaxy 
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Understanding galaxy mergers and their connection 
to these observables requires comparison with theoreti- 
cal predictions. In simulations, galaxies are identified with 
subhalos, the sub s tructures of dark matter halos (e.g . 
iGhigna et all Il998l : iMoore et all Il999l : iKlvpin et all 1 19991 : 



Henceforth, we will use the term galaxy and subhalo 
interchangeably. 

A subhalo forms when two halos collide and a rem- 
nant of the smaller halo persists within the larger final halo. 
Thus, subhalo merger rates are sometimes inferred from halo 
merger ratefl or subhalo distributions, using a dynamical 



^ A subhalo can comprise an entire halo if there are no other 
subhalos within the halo, see [[2] 

^ Halo merger counts and rates have been studied in a vast lit- 
erature, both estimated analytically (e.g., iKaufFmann fc White! 
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friction model to estimate the infall time of satellite galax- 
ies to their halo' s central galaxy (sev eral of these methods 
are compared in i Hopkins et al] |200^. However, a detailed 
understanding of galaxy mergers requires a sufficiently high- 
resolution simulation that can track the evolution and coa- 
lescence of subhalos directly. 

Here we use high resolution dark matter simulations 
to examine subhalo merger rates, counts, and types, their 
mass and redshift dependence, and their relation to their 
host halos. Under the assumption that galaxies populate the 
centres of dark matter subhalo potential wells, our subhalos 
are expected to harbor massive galaxies (L > L*). Although 
our subhalo mass assignment is motivated by semi-analytic 
arguments, our results are independent of any specific semi- 
analytic modeling prescription. 

Previous work on subhalo me rgers include s both 
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_Simha ot al. 2008). Many of these earlier studies concern 
subhalo mergers within a single object (such as the Milky 
Way or a galaxy cluster), others use lower resolution. We 
look at relatively large simulation volumes (100 and 250 

Mpc boxes) with high spatial and temporal resolution. 
In addition, many previous works focus on subhalo mass 
loss and survival rate, while our main interest here is the 
population of resu lting merged subha l os itself (most similar 
to the works of iMaller et al.1 I2OO6I: IGuo fc White! l2008l : 
ISimha et all I2OO8I : lAngulo et aLllioOSh . We characterize 
subhalo merger properties in detail, including satellite merg- 
ers, at high redshift (z = 0.6 — 5), during the peak of merger 
activity. We also investigate the satellite halo occupation 
(number of satellites per halo), and its evolution as shaped 
by the relative merger rates of subhalos and halos. The 
satellite halo occupation is a key elem e nt in the halo model 
llSeliakI [2OO0I: iPeacock fc SraithI I2OO0I : iBerhnd fc Weinberg! 
!2002l : !Coorav fc Sheth|[20o3 ). a framework which describes 
large-scale structure in terms of host dark matter halos. 

In !j2]we describe the simulations, our definitions of ha- 
los and subhalos, their infall mass, and the evolution of the 
satellite fraction. In SjSlwe define our merger criteria and ex- 
amine subhalo merger rates and their relation to halo merger 
rates. In 331we explore the relative contributions of satellites 
and centrals to the merger population, both in terms of par- 
ents and merged children. iJSl gives the distribution of the 
number of mergers within a given look-back time and the 
fraction of halos that host subhalo mergers. In !j6l we show 
how the difi'erence in halo and subhalo merger rates con- 
tributes to the evolution of the satellite halo occupation. 
We summarize and discuss in SjTj The Appendices compare 
subhalo infall mass to infall maximum circular velocity and 
its evolution with redshift, and give our satellite subhalo 
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mass function. The halo occupation and clustering of h igh 
redshift galaxy mergers is treated in !Wetzel et all (|2009! ). 



2 NUMERICAL TECHNIQUES 
2.1 Simulations 

We use two dark matter only N-body simulations of 800^ 
and 1024^ particles in a periodic cube with side lengths 
lOO/i^^ Mpc and 250h~^ Mpc, respectively. For our ACDM 
cosmology (fi^ = 0.25, = 0.75, h = 0.72, n = 0.97 and 
gg = 0.8), in agreement with a wide arra y of observations 
JSmoot et al.l ! 19921: FTegrnark et al.! !2006! : !Reichardt et al.! 
!2008l : !Komatsu et ahlboogl ). this results in particle masses of 
1.4 X lO^h-^MQ (1.1 X IO^H-'^Mq) and a Plummer equiv- 
alent smoothing of 4/i~^ kpc {9h~^ kpc) for the smaller 
(larger) simulation. The initial conditions were generated 
at z = 200 using the Zel'dovich approximation applied to 
a regular Cartesian grid of particles and then evolved using 
the TreePM code described in White (200^) (f or a compari- 



son w ith other codes see !Heitmann et al.ll2008! : !Evrard et al.! 
!2008h . Outputs were spaced every 50 Myr (~ 100 Myr) for 
the smaller (larger) simulation, from 2: '--^ 5 to 2.5. Addi- 
tional outputs from the smaller simulation were retained at 
lower redshift, spaced every ~ 200 Myr down to z = 0.6, 
below which we no longer fairly sample a cosmological vol- 
ume. For mergers, we restrict these later outputs to z < 1.6 
based on convergence tests of the merger rates (we lack suf- 
ficient output time resolution in the intervening redshifts to 
properly catch all mergers; see end of i]3.1l for more). Our 
redshift range oi z = 0.6 — 5 allows us to examine subhalos 
across 7 Gyr of evolution. 

To find the subhalos from the phase space data we 
first generate a ca talog of halos usin g the Friends-of- Friends 
(FoF) algorithm jPavis et al.!! 19851 ) with a linking length 
of 6 = 0.168 times the mean inter-particle spacing. This 
partitions the particles into equivalence classes by linking 
together all particles separated by less than b, with a den- 
sity of roughly p > 3/(27r6'^) ~ 100 times the background 
density. The longer linking length of & = 0.2 is often used. 
However, this linking length is more susceptible to joining 
together distinct, unbound structures and assigning a halo 
that transiently passes by another as a subhalo. Thus, we 
use a more conservative linking length, which for a given 
halo at our mass and redshift regime yields a ~ 15% lower 
mass than b — 0.20 We keep all FoF groups with more than 
32 particles, and we refer to these groups as "(host) halos". 
Halo masses quoted below are these FoF masses. 

When two halos merge, the smaller halo can retain 
its identity as a "subhalo" inside the larger host halo. We 
identify subhalos (and sometimes subhalos within subha- 
los) using a new implementation of the SubGnd algorithm 

^ Many Millennium subhalo studies use Spherical Overdensity 
(SO) halos based on an FoF(0.2) catalogue, in part to take out 
the ex tra structure joined by the larger linking length. [Neto et al] 
1I2OO7I') compares the FoF(0.2) halo centres and those for the SO 
halos use d to define the corresponding subhalo populations. More 
generall y. IWhitd (I2OOII') compares different mass definitions in de- 
tail, and !Cohn fc Whitel bOOSl ) discusses the rel ation of FoF(0.2) 
FoF(0.168), SO(180) and Sheth- Tormen ( Shcth fc Tormenlll999l 
based on 6 = 0.2) masses at high redshift. 
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l|Springel et al.ll200ll ). We take subhalos to be gravitation- 
ally self-bound aggregations of particles bounded by a den- 
sity saddle point. After experimentation with different tech- 
niques we find this method gives a good match to what 
would be selected "by eye" as subhalos. We use a spline 
kernel with 16 neighbours to estimate the density and keep 
all subhalos with more than 20 particles. The subhalo that 
contains the most mass in the halo is defined as the central 
subhalo, all other subhalos in the same halo are satellitesQ 
The central subhalo is also assigned all matter within the 
halo not assigned to the satellite subhalos. The position of 
a subhalo is given by the location of its most bound parti- 
cle, and the centre of a host halo is defined by the centre 
of its central subhalo. For each subhalo we store a number 
of additional properties including the bound mass, velocity 
dispersion, peak circular velocity, total potential energy, and 
velocity. 

Figure [T] (top) shows the projected image of a sample 
halo of mass 2.2 x W^^ Mq, which hosts 12 satellite sub- 
halos, at z — 2.6 in the 100/i~^ Mpc simulation. (At this 
redshift in our simulations there are approximately 200 ha- 
los per (100/i~^ Mpc)^ above this mass.) Halos at this mass 
and redshift regime are dynamically active and often highly 
aspherical, with mean axial ratio (ratio of smallest-to-largest 
semi-major axes) of ~ 0.5, and r ecently merged halos have 
even more discrepant axial ratios (jAllgood et al.|[200d ) . Fig- 
ure [T] also shows how a halo's densest region (figure centre), 
centre of mass, and substructure distribution can all be off- 
set from one another. Because of these asymmetries, satellite 
subhalos sometimes extend well beyond r2ooc, which we will 
call the halo virial radiuslf] 

2.2 Subhalo Tracking 

We identify, for each subhalo, a unique "child" at a 
later t i me, using subhalo tr a cking si milar to Springel et al.l 



20051). iFaltenbacher et all (|2005| ). lAllgoodI (|2005l ') and 



iHarker et al. I l|2006l ). We detail our method to illustrate the 



subtleties which arise and to allow comparison with other 
work. 

We track histories over four consecutive simulation out- 
puts at a time because nearby subhalos can be difficult to 
distinguish and can "disappear" for a few outputs until their 
orbits separate them again. For each subhalo with mass Mi 
at scale factor ai, its child subhalo at a later time (scale 
factor a2 and mass A/2) is that which maximizes 



a = /(Mi,M2)ln-M ^K]<^?, 



where 



/(Ml, Ma) 



1 - \Ml-M2\ ^ 

1 Ml > M2 



(1) 



(2) 



In most cases, the central subhalo in SubBnd is built around 
the most bound and most dense particle in the group. However, 
this is not always the case. 

^ The halo virial radius, raooci i-e- the radius within which the 
average density is 200 X the critical density, is calculated from 
the FoF (b = 0.16 8) mass by first con verting to M20OC assuming 
a spherical NFW llNavarro et aLlFiggQ) density profile, and then 
taking Afaooc = 200^pc4„o^. 
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Figure 1. Top: Projected image of a halo of mass 2.2 X 
10^^/i~^A/q at 2 = 2.6 which hosts 12 satellite subhalos. Parti- 
cles assigned to the central subhalo are red, while those assigned 
to satellite subhalos are blue. Dot-dashed circle shows the halo's 
virial radius (r200c), derived from its mass assuming a spherical 
NFW density profile, while the solid circles highlight the satellite 
subhalos and scale in radius with their mass. The central sub- 
halo has bound mass of 1.8 X lO^^h-^MQ, so nearly 20% of the 
halo's mass lies in satellite subhalos. Bottom: Tracking histories 
of massive satellite subhalos in the above halo. Large dots show 
the positions of satellite subhalo centres at z = 2.6 for subha- 
los that had a mass > l0^^h~^MQ when they fell into the halo. 
Small dots show their positions (relative to that of the halo cen- 
tre) at each output (spaced 50 Myr) back 800 Myr. Thin blue 
curves show subhalo trajectories when they are satellites while 
thick red curves show when they are centrals (before falling into 
the halo). Dotted lines indicate when the parent-child assignment 
has skipped an output (during a fly- by near another subhalo). 
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and where (pu is the potential of particle i computed using 
all of the particles in subhalo 1, and the sum is over those of 
the 20 most bound particles in the progenitor that also lie in 
the candidate child. We track using only the 20 most bound 
particles since our ultimate interest is in galaxies, which we 
expect to reside in the highly bound, central region of the 
subhalo (20 is the minimum particle count for our subhalos). 
We do not use all the progenitor particles because summing 
over all of the particles in the progenitor that also lie in the 
child candidate leads to instances of parent-child assignment 
in which the child subhalo does not contain the most bound 
particles of its parent. Finally, we weight against large mass 
gains with a mass weighting factor so that smaller subhalos 
passing through larger ones and emerging later on the other 
side are correctly assigned as fly-by's and not mergers. We 
find ~ 95% of subhalos have a child in the next time step, 
with ~ 4% percent skipping one or more output times and 
~ 1% having no identifiable child. 

Figure [T] (bottom) shows the tracking histories of the 
most massive subhalos within the halo shown in Fig. [1] (top), 
for subhalos that had a mass > W^^ Mq when they fell 
into the halo. The upper-most subhalo was a separate halo 
(M = 1.3 X lO^^/i"^M0) hosting its own massive satellite 
subhalo. It then fell into the main halo and then both (now 
satellite) subhalos merged with each other. The right-most 
subhalo is an example of two separate halos falling into the 
main halo, becoming satellites, and subsequently merging 
with each other. Finally, the track through the centre shows 
a single subhalo falling towards the central subhalo. Instead 
of merging with the central, it passes through as a fly-by. 

All of these tracks show instances where the parent- 
child assignment algorithm has skipped an output (dotted 
lines) as one subhalo passes through another and re-emerges 
on the other side. Note also that, while the dot-dashed circle 
shows the halo virial radius, r2ooc (at the last output), based 
on its mass assuming a spherical NFW density profile, the 
locations of the transitions of central subhalos to satellites 
during infall show that the spherical virial radius is only a 
rough approximation. 

As a halo falls into a larger host halo and becomes 
a satellite subhalo, mass loss from tidal stripping can be 
extreme (90% or more) as the satellite subhalo orbits to- 
wards the centre of its host halo. Figure [2] shows the evo- 
lution from infall of a long-lived satellite with infall mass 
Mini = 6 X 10^^ Mq and maximum circular velocity at 
infall of Vc.inf = 258 km/s. Mass and maximum circular 
velocity (K"!max) are stripped by ~ 90% by the time the 
satellite first reaches pericentre. After pericentric passage, a 
satellite can also gain mass and circular velocity as it moves 
away from the dense centre of its host halo. This is shown 
in Fig. [2] where mass and circular velocity strongly corre- 
late with radial distance throughout the orbit. This is par- 
tially an effect of our subhalo finder: the density contrast 
of a subhalo, which defines its physical extent and hence its 
mass, drops as it moves to a region of higher background 
densitylfl However, this is also driven by a physical effect: 
as a satellite approaches the dense centre of its host halo, 
it will be compressed by tidal shocks in response to the 
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Figure 2. Top: Evolution of mass (solid curve) and circular 
velocity (dashed curve) as a function of time since infall for a 
satellite that fell into a halo of mass 8 X 10^^h~^ Mq at ^ = 3.7. 
Bottom: Radial distance of satellite from centre of host halo. 
Mass and circular velocity exhibit correlations with radial dis- 
tance, such as mass gain as the satellite recedes from the centre 
of its host halo. This satellite experienced no major merger activ- 
ity throughout its history. 



rapidly increasing potential (see Fig. 12 of Diemand et al.l 
l2007l : iGnedin fc Ostrikedl 19971 : iDekel et al.ll2003l ). 

Since a central subhalo is defined as the most massive 
subhalo, a satellite subhalo can become a central (and vice 
versa) , which we refer to as a "switch" We find that ~ 4% 
of all centrals at any output become satellites within the 
same halo in the next output, and of these a third im- 
mediately return to being centrals in the following output. 
Switches are twice as common during a merger event, ei- 
ther when two satellites merge, or when a satellite merges 
with the central, forming a less dense object and allowing 
another satellite to bind more mass to itself. When a satel- 
lite switches to a central and back to a satellite, the original 
satellite can be mistakenly assigned as a direct parent of the 
final satellite (and thus the central is assigned no parent) 
since our child assignment weights against large mass gains 
(which occurs when a satellite becomes a central). We fix 
these distinct cases by hand. 

These switches highlight the fact that the distinction 
between a central and satellite subhalo is often not clear- 
cut: at this redshift and mass regime massive halos undergo 
rapid merger activity and thus often are highly disturbed 
and aspherical, with no well-defined single peak that repre- 
sents the center of the halo profile. 



In more detail, a switch occurs when the density peak of a satel- 
^ A fly-by, i.e. when a subhalo passes through and is temporarily lite (above the background) contains more mass than is within the 

indistinguishable from a larger subhalo, is an extreme case of this. central subhalo's radius at the position of the satellite 
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2.3 Subhalo Mass Assignment 

Since we use subiialos as proxies for galaxies, we track sub- 
halo mass that is expected to correlate with galaxy stellar 
mass. Galaxies form at halo centres as baryons cool and adi- 
abatically contract toward the minimum of the halo's poten- 
tial well, which leads to a correlation between halo mass and 
galaxy stellar mass llWhite fc Reeslll978l : iBlumenthal et alj 
ll986l : lDubin"skilll994l : lMo et aLlllQQsl ). When a halo falls into 
a larger halo and becomes a satellite subhalo, its outskirts 
are severely stripped as discussed above, but its galaxy's 
stellar mass would be little influenced as the galactic ra- 
dius is typically ~ 10% that of the subhalo radius. This 
motivates assigning to subhalos their mass at infall, Minf, 
which is expected to correlate with galactic stellar mass 
throughout the subhalo's lifetime. The subhalo infall mass 
function has been successful at reproducing the observed 
galaxy luminosity function and clusterin g at low redshift s 
(jVale fc Ostrikerli2006l : rWang et al. 2006: Yang et al.ll2009h . 
Maximum circular velocity at infall, K.in f , has also been suc- 
cessfully matched t o some observations ( Conrov et al.|[2006l : 
iBerrier et al.llioO^ '). In Appendix A we show the relation 
between Mini and Vc.inf and its redshift evolution. 

Our prescription for assigning Minf to subhalos is as fol- 
lows. When a halo falls into another and its central subhalo 
becomes a satellite subhalo, the satellite is assigned Mini 
as the subhalo mass of its (central) parentlf] If a satellite 
merges with another satellite, the resultant child subhalo 
is assigned the sum of its parents' Minf. Since the central 
subhalo contains the densest region of a halo, inter-halo gas 
is expected to accrete onto it, so we define Mi^i for a cen- 
tral subhalo as its current self-bound subhalo mass, which 
is typically ~ 90% of its host halo's mass|f| However, since 
a central subhalo can switch to being a satellite, while a 
satellite switches to being a central (all within a single host 
halo), we require an additional rule because using the above 
simple assignment of Minf to centrals would lead to a cen- 
tral and satellite in a halo each having the halo's current 
bound mass. Thus, we assign a central to have Minf as its 
current self-bound subhalo mass only if it was the central 
in the same halo in the previous output. Thus, if a satellite 
switches to a central and remains the central for multiple 
outputs, it has robustly established itself as the central sub- 
halo, so it is assigned its current self-bound mass. However, 
if a central was a satellite (or a central in another smaller 
halo) in the previous output, it is assigned the sum of its 
parents' Minf, with the additional requirement that its i\finf 
cannot exceed its current self-bound mass. 

A small fraction (~ 4%) of satellites composed of at 
least 50 particles are not easily identifiable with any progeni- 
tor subhalos and thus cannot be tracked to infall. On inspec- 
tion, we find that these "orphaned" satellites are loosely self- 



° There is some dependence on output time spacing in this def- 
inition: since a central subhalo typically continues to gain mass 
before it falls into a larger halo, shorter time steps (which catch 
it closer to infall) lead to a higher subsequent satellite Afinf . Dou- 
bling the output spacing leads to satellites with ^ 10% lower 
infall mass. 

^ Though because of finite resolution, we find a weak systematic 
drop in Mcon/Mhaio with halo mass, varying from 93% to 87% 
for 10^"*^ to lQ^'^h~^ AtQ in the higher resolution simulation. 



bound portions of a central subhalo, remnants from a colli- 
sion between a satellite and its central subhalo that soon re- 
merge with the central. Given their origins and fates, these 
orphans are not expected to host galaxies and are ignored. 

Although we track subhalos down to 20 particles, we 
impose a much larger minimum infall mass to our sample to 
avoid selecting subhalos that artificially dissolve and merge 
with the central too early. This requires sufficient resolution 
of the radial density profile of a satellite subhalo at infall: 
if the satellite's core is smaller than a few times the force 
softening length, its profile will be artificially shallow and it 
will be stripped and disrupted prematurely (see Appendices 
for details). For calibration, we use the regime of overlap 
in mass between our two simulations of different mass res- 
olution, requiring consistent subhalo mass functions, halo 
occupation distributions, and merger statistics for a mini- 
mum infall mass. For example, going too low in mass for the 
larger simulation resulted in more mergers and fewer satel- 
lites than for the same mass range in the smaller simulation. 
In the larger simulation, our consistency requirements led us 
to impose Minf > W^^h~^ Mq. For a fixed number of parti- 
cles per subhalo, this scales down to Minf > 10^^ Mq in 
the smaller, higher resolution simulation. Note that halos of 
mass 10^^ {10^^)h~^ Mq cross below Af, , the characteristic 
mass of collapse, at z = 1.5 (2 = 0.8), so we probe massive 
subhalos across most of our redshift range. 

At 2 = 2.6, there are ~ 16000 subhalos with Minf > 
lO"/i"^M0 in our lOO/i"^ Mpc simulation and 9400 sub- 
halos with Minf > 10^^h~^MQ in our 250/i~^ Mpc simu- 
lation. At 2 = 1, there are ~ 29400 (2500) subhalos with 
Minf > 10" (10")/i"^MQ in our lOO/i"^ Mpc simulation. 



2.4 A Note on Stellar Mass and Gas Content of 
Subhalo Galaxies 

A galaxy's stellar mass is expected to be a non-linear, 
redshift-dependent function of its subhalo mass. An approx- 
imate relation based on abundance matching is given in 
IConrov fc Wechsleij (120081 ). At 2 = 1, they find that subha- 
los of infall mass 10" (10^^)ft"^MQ host galaxies of stellar 
mass ~ 10^ (10^° '^)Mq. At 2 = 2.5, subhalos at the above 
masses are expected to host lower mass galaxies, though 
quantitative relations at this redshift are less certain, and 
our subhalo and halo finders differ in detail from theirs. Our 
sample of Afinf > 10^^ Mq subhalos approximately cor- 
responds to L > L, galaxies at the redshifts we examine. 

Although our simulations do not track the baryonic con- 
tent of subhalos, most massive galaxies are gas-rich at high 
redshift. That is, at 2 » 1, 70% to 90% of L ~ Lt galaxies ar e 
observed to be blue jCooper et al.ll2007l : [Cerke et al.ll2007l '). 
possessing enough gas to be actively star forming. The frac- 
tion of gas-rich galaxies at higher redshift i s more poorly 
const rained but is thought to be higher (|Hopkins et al.l 
2008). Thus, we anticipate that most, if not all, mergers 
we track have the capacity to drive galaxy activity such as 
starbursts and quasars. 



2.5 Satellite Fraction 

To frame our ensuing discussion of mergers, and the relative 
importance of satellite and central subhalos. Fig. |3] shows 
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Figure 3. Evolution of tiie fraction of all subhalos that are 
satellites for Afj^f = 10^^ - 10^^h~^MQ (triangles) and Mj^f = 
10^2 - lO"/i-iM0 (squares). Das hed lines show fits to Eg. [H 
Dotted line shows fit to Eq.|4](that of lConrov fc Wechslej ll2008l )'). 



the evolution of the satellite fraction, nsatciiito/fisubhaio, for 
subhalos of a fixed masj^. which grows monotonically with 
time from z — 5 — 0.6, as 



^satellite 



= C-e-— (3) 

' ''subhalo 

with /3 — 9.7 valid across all subhalo masses and C — 
0.26 (0.24) for lower (higher) mass. The satellite fraction 
decreases with increasing subhalo mass since more massive 
subhalos are more likely to be centrals, and for a fixed sub- 
halo mass the satellite fraction increases with time as the 
number of high mass halos hosting massive satellites in- 
creases. The increase in the satellite fraction is slowed at 
late times because the number of massive satellites in halos 
of a fixed mass decreases with time as the satellites coalesce 
with the central subhalo (see |6]). 

Note th at Eq.[3]predicts a m u ch hig her satellite fraction 
than that of lConrov fc Wechslerl |200i), who found 

Jhs^ = 0.2 - —z. (4) 

^subhalo 3 

It is unlikely that the difference is driven by numerical ef- 
fects, since both their and our simulations are of similar mass 
resolution and volume, and both analyses are based on simi- 
lar subhalo infall mass cuts (we see similar satellite fractions 
selecting instead on fixed Vcinf)- One likely factor is that, 
as noted above, our smaller output spacing yields higher in- 
fall mass for satellite subhalos, since we catch halos closer 
to infall when their mass is higher. In additio n, they use a 
different halo and subhalo finding algorithm (|Klvpin et al.l 
ll999l : lKravtsov et al.ll2004l '). and their higher as (0.9 rather 
than our 0.8) would give a lower merger rate and thus fewer 
satellites (one expects satellite survival timescales not to 
change). If normalized to the same satellite fraction at a 
given epoch, the redshift evolution of Eqs. O and [3] are in 
rough agreement. 

We show the satellite fraction for two mass bins, but since the 
mass function falls exponentially at these masses and redshifts, 
almost all objects are at the low end of the mass bin. Using instead 
a minimum mass cut changes our results by only a few percent. 



3 MERGER CRITERIA AND RATES 
3.1 Merger Criteria 

For two parents with A/inf,2 < A/inf,i sharing the same child 
subhalo at the next output they appear, a child is fiagged 
as a (major) merger if ^^j'"^;^ > i. If a child has more than 
two parents, we count multiple mergers if any other parents 
also exceed the above mass ratio with respect to the most 
massive parent. 

Our mass ratio represents a trade-off between strong 
mergers (to maximize signal) and frequent merging (for 
statistical power). Galaxy mergers with stellar mass ratios 
closer than 3:1 are expected to drive interesting activity, e.g. 
quasars and starbursts, as mentioned abover^l 

For generality, the distribution of merger (infall) mass 
ratios, R:l, for both halos and subhalos in the mass and 
redshift regimes we consider can be approximated by 



fiR) (X R- 



(5) 



in reasonable agreement with the R ^' distribution of 
galaxy mass ratios at 2 < 0.5 f ound in hydrodynamic sim- 
ulations by iMaller et al] l|2006l ')Fl Thus, the counts/rates 
of mergers with (infall) mass ratio closer than R:l can be 
approximately scaled from our results through the relation 



iV(< R) = 8.6 {R° ' - l)iV(< 3). 



(6) 



Sometimes me rger criteria ar e bas ed on instantaneous 
subhalo mass gain iThacker et "all (e.g., [20061 ). However, as 
exemplified in Fig. [2j subhalos can gain significant mass 
without coalescence. We find that most cases of significant 
subhalo mass gain are not a two-bod y coalescence, and so we 
do not use this to select mergers f seelWetzel et al.l (1200 9) for 
more detail, also related results in lMaulbetsch et al.l "i 2007i )). 

Sufficient output time resolution is necessary to prop- 
erly resolve the merger population. Unless otherwise stated, 
we use the shortest simulation output spacing to define 
the subhalo merger time interval, corresponding to 50 Myr 
(~ 100 Myr) for Mi„f > lO" {\{i^'^)h-^ Mq at 2 > 2.5, 
and ~ 200 Myr for all masses at z < 1.6. By tracking par- 
ent/child subhalo assignments across multiple output spac- 
ings, we found that longer output time spacings result in 
a lower subhalo merger rate, arising from a combination 
of effects. As stated above, shorter time steps catch halos 
closer to their infall, resulting in satellites with higher Minf, 
and hence more major satellite-central mergers. In addition, 
satellite-satellite mergers can be missed if their child merges 
into the central subhalo before the next time step. How- 
ever, smaller output spacings are also more susceptible to 
switches, which can yield satellites with higher mass. 

For halos, very short time steps can artificially enhance 
the merger rate, catching transient behavior as halos inter- 
sect, pass through each other, and re-merge. While we use 



Since galaxy stellar mass is a non-linear function of subhalo 
mass, a subhalo mass ratio of 3:1 may correspond to a galaxy 
merger of a more or less discrepant mass ratio. However, since the 
A^stcllar - A^s ubhalo relation is expected to peak for Mgubhalo ~ 
IQI^/i-IMq l lConrov fc Wechsle Jr200^ . we do not expect this 



effect to strongly bias our results. 

^■^ The distribution of merger mass ratios also agrees well with 
the fit for halos at 



I merger mass ratios also ag r ees w e, 
: provided bv lWetzel et al'l ||200^ . 
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FoF(0.168) instead of FoF(0.2) halos to help minimize ar- 
tificial bridging effects, for a 50 Myr time spacing halo re- 
mergers still constitute a significant fraction of halo mergers. 
Thus we calculate all our halo mergers using > f 00 Myr out- 
put spacings, for which halo re-mergers constitute only a few 
percent of all mergers. 

Given our simulation volumes we are unable to probe 
statistically significant subhalo merger counts in halos more 
massive than 4 x 10^^/i"^Mq at z > 2.5 and M ~ 
10^"^ A4q aX z < 1.6. Ifowever, the majority of recently 
merged subhalos are centrals (discussed below), and select- 
ing on higher mass host halos for a fixed subhalo mass re- 
stricts increasingly to recently merged satellites (of which 
there are comparatively fewer) . Furthermore, nearly all sub- 
halo mergers at a given mass occur in halos at most 5 — lOx 
times more massive. 



3.2 Fits to Simulation 

Figure |4] shows the evolution of the merger rate per object — 
defined as the number of mergers per time per object — for 
both subhalos and halos of the same (infall) mass|^ We fit 
the merger rate per object as 



A{1- 



(7) 



where rimcrgc is the number of (sub)halos whose parents 
match our mass ratio selection, Jiobj is the total number 
of (sub) halos within the same mass range at the same out- 
put, and dt it the time interval between consecutive outputs. 
The best-fit values in each redshift regime are shown in Ta- 
ble[T] Note that the merger rate we examine is the number of 
mergers per time per object, different from another common 
definition, the number of mergers per time per volum^^ 

The relation of subhalo mergers to halo mergers is non- 
trivial. Although subhalo mergers are the eventual result of 
halo mergers, the former are governed by dynamics within 
a halo and the latter by large-scale gravitational fields. For 
halos, both the slope and the amplitude of the merger rate 
exhibit little dependen ce on mass or redshift. Similarly, e.g. 
iFakhouri fcl^l|2008l ) found a = 2 - 2.3 for all halo masses 
at 2 < 6, and weak halo mass dependence of the amplitude. 
In contrast, the subhalo merger rate amplitude has strong 
dependence on mass and its slope depends strongly on red- 
shift. Relative to the halo merger rate, the subhalo merger 
rate is lower in amplitude than that of halos of the same (in- 
fall) mass, and, most notably at z > 2.5, the rate of subhalo 
mergers falls off significantly more sl owly that that o f halos. 
This is consistent with earlier work: iDe Lucia et al.l (|2004| ) 



found a higher merger fraction for halos than subhalos, and 
IGuo fc Whit3 (2008) found strong mass dependence of the 
galaxy merger rate amplitude and that the slope becomes 
much shallower at z > 2 (see also iMateu s 200.8). (Note that 
as our halos are FoF(0.168) halos and so our merger rates can 
differ from those for FoF(0.2) halos; merging occurs sooner 
for a finder with a larger linking length.) 

We now focus on the relation between halo and subhalo 
merger rates to understand these trends with time. 

3.3 Subhalo vs. Halo Merger Rates 

A simple analytic argument based on dynamical infall time, 
i.e., subhalo mergers are simply a delayed version halo merg- 
ers, leads one to expect that subhalo merger rates simply 
track those of halos: they have the same time evolution, with 
the subhalo merger rate having a higher amplitude. In this 
argument, when two halos merge, the new satellite galaxy 
collides with the other central galaxy within a dynamical 
friction timescale, approximated by 

Mhalo/Msat 



where tdyn 



*'""^^^^°in(l + A4a,o/M,at) 
O.liHubble, tHubble 



tdy 



(8) 



accounts f or the ensemble av e raged satellite orbital pa- 
rameters (jConrov e t aP '200 71: iBinnev fc Tremaind l2008l : 
iBovlan-Kolchin et al. 2008 : Jiang et al.l 120081 ). Thus, for 
a fixed mass ratio, letting nio = Mhaio/Afsat/ ln(l -I- 

Mhalo/Msat), 



sat, merge 



0.1 Comotuub 



0.1 Comot. 



(9) 



The evolution of the halo merger rate per object during 
matter-domination (valid at the high redshifts we examine), 
where a <x t^^^ , is approximately 



nhaloCtt \ t 



(10) 



with t-M some proportionality constant. Assuming all halo 
mergers lead to satellite-central subhalo mergers on a dy- 
namical friction timescale, the subhalo merger rate per ob- 
ject would evolve with time as 



^subhalo 



= (1- 0.1 Como)-3" ( _ 



(11) 
(12) 



and so the subhalo merger rate would simply track that of 
halos of the same mass, but with higher amplitude. 



Merger fraction errors are calculated using binomial statistics. 
Given M mergers out of objects, the mostly likely fraction is 
t n^lM -iU ■ 2 M(N-M) + 1+N 

f = M/N with variance aj = (jv+2)S(iV+3) ' 

The former is simply the latter divided by the number density 
of objects of the same mass, but the merger rate per volume has 
qualitatively different behaviour because the (comoving) number 
density of objects at a fixed mass increases with time in a redshift- 
dependent manner. Specifically, at high redshift where the mass 
function rapidly increases, the merger rate per volume for both 
halos and subhalos increases with time, reaching a peak at 2 ~ 
2.5. Below this redshift, it decreases with time in a power law 
manner as in Fig.|4l though with a shallower slope of q 1.5. 



3.4 Resolving the Discrepancy 

As Fig. |4]shows, however, this tracking does not occur, par- 
ticularly at high redshift where the slope of the subhalo 
merger rates is much shallower than that of halos. One rea- 
son for this is that we compute the merger rate per object, 
in which we divide by the number of objects at the given 
mass, riohiim). Since halo masses are added instantaneously 
during halo mergers, a recently merged halo will instantly 
jump to a higher mass regime (with smaller riobj in Eq. [7}. 
In contrast, the central subhalo of the resultant halo will 
remain at a smaller mass for some time until its mass grows 
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Figure 4. Left: Merger rate per object of subhalos (solid points) and halos (open points) as a function of the scale factor, for 
Minf = 10" -lOl2/i-lM0 (triangles) and Mj^f = lO^^ - IqIS/i-IMq (squares) at 2 = 5-2.5. Solid (dashed) curves show fits of subhalo 
(halo) merger rates to Eq. [7] with fit parameters given in Table[T] Right: Same, but at z = 1.6 — 0.6. Also shown at left is the fit to the 
halo merger rate for halos with a central subhalo of Mjnf = lO" — lQ^^h~^ Mq (dotted). 
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10" - 10^2 


10^2 - 10^3 
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10^2 - 10^3 


Halos 
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0.029 


0.032 


0.034 


0.034 












a 


2.3 


2.3 


2.0 


2.2 


Subhalos 


A 


0.093 


0.065 


0.016 


0.016 




a 


1.1 


1.1 


2.1 


2.6 



Table 1. The amplitude. A, and power law index, a, for halo and subhalo merger rates fit to Eq. [7] 



from stripping of the new satellite subhalo (i.e., recently 
merged halos have lower than average Minf^ccn/Afhaio), so 
the central's riobj is higher. 

In Fig. |4] (left) the dotted black curve fits the halo 
merger rate for halos selected with the same mass cut on 
their central subhalo infall mass. The amplitude is signifi- 
cantly smaller, and the slope is also shallower (a — 1.6) than 
for halos selected on full halo mass, showing that the above 
effect is stronger at earher times. By the argument above, 
the ratio of amplitudes of the merger rate per object of ha- 
los to subhalos is Whaio/'^cen. A fixed (sub)halo mass cut 
probes lower M/M,{t) at later times, and the mass func- 
tion drops exponentially with increasing M/M,{t) at these 
masses. Thus, a fixed A/inf,con/Afhaio after a halo merger 
means nhaio/jicon becomes closer to unity at later times, 
leading to the shallower slope. 

The measured subhalo mergers have an even lower am- 
plitude and shallower slope than the dotted line in Fig. U 
driven by two additional effects. First, a halo major merger 
might not lead to a subhalo major merger since the satellite- 
central merger mass ratio can be smaller than the mass 
ratios of their source halos. This is because Mint natu- 
rally grows for a central but only grows for a satellite 
if it has a merger before coalescing with the central (see 
I Wang fc Kauffmannlliooi . for a detailed analysis of this ef- 
fect in terms of assigning baryons to subhalos). Since halos 
grow in mass more quickly at higher redshift, this effect 
is stronger at earlier times, further flattening the subhalo 



merger rate slope. Second, there is a significant contribu- 
tion of recently merged satellites to the merger population. 
We find that satellites are twice as likely to have had a re- 
cent merger as centrals of the same mass, regardless of mass 
cut and redshift (see 311 for more detail). This enhances the 
merger rate, with a stronger enhancement at later times 
since the satellite fraction grows with time as in Fig. |3l 

At lower redshift. Fig. |4] (right) shows that the ampli- 
tude of the subhalo merger rate remains lower than that 
of halos, but the slopes become similar, indicating the ef- 
fects examined above become less time-dependent. Since the 
masses we probe are crossing M*(t) as these redshifts, the 
reduced amplitude from subhalo vs. halo mass cut becomes 
less sensitive to time. Similarly, our mass range crossing 
Mt{t) means that halo mass growth slows, so the fraction 
of halo major mergers that leads to subhalo major merg- 
ers remains roughly constant with time. Finally, as shown 
in Fig. [S] the satellite fraction growth asymptotes at lower 
redshift, which means that the enhancement from recently 
merged satellites remains roughly constant. 



4 SATELLITE VS. CENTRAL MERGERS 

The stereotypical galaxy merger is a satellite coalescing with 
the central in its halo and producing a central merger rem- 
nant. These mergers do dominate the merger population, 
both in parent types (central-satellite) and child type (cen- 
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Figure 5. Satellite merger fraction at 2: = 5 — 2.5 (left) and 2 = 1.6 — 0.6 (right) for Mj^f = 10^^ — 10^^h~^ Mq (triangles) and 
Afi„f = 10^^ — 10^^ Mq (squares). Closed point/solid lines indicate the fraction of all subhalo mergers that result in a satellite, while 
open points/dashed lines indicated the fraction of all subhalo mergers that arise from satellite-satellite parents. 



tral). However, while satellites form the minority of the sub- 
halo population at all epochs (see Fig.[5J, satellites are twice 
as likely to have had a recent merger as centrals of the same 
mass, regardless of mass and redshift. 

Since the identities of satellite vs. central subhalos 
at this mass and redshift regime are not clear-cut (from 
switches), we characterize mergers both in terms of their 
parent types (central/satellite) and resulting child types. For 
mergers resulting in centrals, this ambiguity is not impor- 
tant: 97% arise from satellite-central parents, while the other 
3% arise from satellite-satellite parents during switches 

Recently merged satellites are a more varied population. 
Figure [5] shows the contributions of satellite mergers to the 
overall merger populations as a function of the scale factorF^ 
The fraction of all subhalo mergers that result in a satellite is 
- 30% for Minf = 10" - 10^^/i"^MQ, with little dependence 
on redshift. For Minf = 10^^ - lO^^/i"^M0 it is 15 - 20% 
at z ^ 2.5 and rises to ~ 35% at 2 < 1.6. Thus, at lower 
redshift {z < 1.6) where the satellite fraction asymptotes 
to ^ 25%, the fraction of mergers that result in a satellite 
roughly reflects the satellite fraction as a whole. 

Of these recently merged satellites, 20% — 35% come 
from satellite-satellite parents within a single halo, while 
~ 7% arise when a central-satellite merger occurs in a halo 
as it falls into a larger halo, becoming a satellite. The rest 
arise from switches, i.e. a satellite merges with a central, and 
the resulting subhalo no longer is the most massive subhalo, 
thus becoming a satellite. These switches occur primarily 
in halos only a few times more massive than the satellite, 
typically for satellites in close proximity to their central. 
At higher halo masses, recently merged satellites are dom- 
inated by satellite-satellite parents. These satellite-satellite 
mergers preferentially occur in the outer regions of a halo 
and are comparatively less common in the central regions. 

A few percent arise from central-central parents, when the 
central regions of two halos coalesce so quickly that they are not 
seen as satellite-central subhalos given finite time resolution. We 
include these as satellite-central parents. 

Satellite merger fractions are boxcar-averaged across three 
consecutive outputs to reduce noise from small number statistics. 



We examine in more detail the en v ironm ental dependence 
of subhalo mergers in IWetzel et al] (|2009l ). 

Considering instead only parent types, Fig.[5]shows that 
5—10% of all mergers come from satellite-satellite parents at 
z ~ 2.5, a fraction which increases to 10 — 15% at z — 1.6 and 
remains flat thereafter. Approximately 80% of all satellite- 
satellite mergers lead to a satellite child, while the rest lead 
to a central during a switch. 



5 GALAXY AND HALO MERGER COUNTS 
5.1 Counts of Recent Mergers 

The distribution of the number of mergers per object within 
a fixed time interval gives the fraction of objects at a given 
epoch that might exhibit merger-related activity or mor- 
phological disturbance 13 Multiple mergers as well might 
contribute to specific properties, e.g ., the formation an d 



mass growth of ellipt ical gala xies (e.g.,|B ovlan-K olchin et al 



20071 ) 



2005 : [Robertson et al. 2006; Naab et al . 2006; C onrov et al 



Figure [6] shows the fraction of subhalos at z — 2.6 and 
2 = 1 with a given number of mergers in the last 1 Gyr. At 
z — 2.6, 30% of subhalos have suffered at least one major 
merger. Interestingly, this fraction is nearly constant across 
the mass regimes we probe. In contrast, the fraction of halos 
with at least one major merger within 1 Gyr is about twice 
as large, with stronger mass dependence: higher mass halos 
experience more mergers. At high redshift, halos are also 
significantly more likely to have undergone multiple mergers 
than subhalos, which builds up the satellite population. 

At 2 = 1, recent mergers of subhalos and halos become 
less common, with only 8% of subhalos having suffered at 
least one major merger in the last 1 Gyr. Objects which have 
had 1 or 2 mergers are still more common for halos than 
subhalos, though high mass subhalos exhibit a much higher 
fraction of 3 or more mergers than halos of the same (infall) 

This differs from the merger rates of ^ since we are tracking 
the histories of individual objects selected at a given redshift. 
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Figure 6. Top; Fraction of subhalos at 2 = 2.6 that have had 
a given number of mergers in the last 1 Gyr, for Mj^f = lO'^^ — 
IQI^/i-IMq (thin lines) and Mi^f = lO^^ _ io13/i-1Mq (thick 
lines). Dashed lines show the same, but for host halos. 30% of 
subhalos have had at least one merger (independent of mass), 
while for halos this fraction is 50% (lower mass) and 57% (higher 
mass). Inset shows detail for > 2 mergers. Bottom: Same, but at 
z = 1. 8% of subhalos have had at least one merger (independent 
of mass), while for halos this fraction is 11% (lower mass) and 17% 
(higher mass). High mass subhalos show a much larger fraction 
of multiple mergers. Inset shows detail for > 1 merger. 



though, as one expects relevant lifetimes to depend strongly 
upon galaxy mass and merger ratio. Observables might de- 
pend upon dynamical time as well, although many quasar 
triggering efltects might be related to microphysics-small 
scale interactions close to the merger-that do not evolve 
with timeH 

Figure [7] shows the evolution of /on for halos hosting 
subhalo mergers within the last 200 Myr. The same quan- 
tity is shown for halos with recent mergers themselves. At 
high redshift, /on for halo mergers shows a steep decline from 
the decreasing halo merger rate. However, /on for subhalo- 
subhalo mergers is fiat from z = 2.5 — 5, because the subhalo 
merger rate per object decreases while the number of mas- 
sive satellites in a given mass halo rises, causing the number 
of massive subhalo mergers within the halo to remain con- 
stant. At low redshift, where the satellite population grows 
more slowly, the evolution of this fraction for subhalo merg- 
ers more closely parallels that of halo mergers. 



6 EVOLUTION OF THE SATELLITE HALO 
OCCUPATION 

The redshift evolution of the satellite galaxy populations 
of dark matter halos is shaped by halo vs. galaxy mergers: 
halo mergers create satellites while galaxy mergers remove 
themB If the infall rate of satellites onto a halo is different 
than the satellite destruction rate, the satellite halo occupa- 
tion will evolve with time. 

As shown in Fig.|4l at z > 2.5, the merger rate of subha- 
los is significantly lower and shallower in slope than that of 
halos, implying that subhalos are being created faster than 
they are destroyed (at a rate decreasing with time). Con- 
versely, a,t z < 1.6 the merger rates of halos and subhalos 
exhibit approximately the same redshift dependence, and 
their amplitudes are similar (also recall from i|3]that not all 
halo major merger lead to subhalo major mergers). Thus, 
for halos of a fixed mass, we expect a rapid rise in the satel- 
lite halo occupation prior to 2 ~ 2 and a levelling-off with 
more gradual evolution at lower redshift. 



mass. (Subhalos that have undergone 1 merger can be either 
satellites or centrals, those that have undergone 2 or more 
mergers are almost entirely centrals.) The build-up of the 
satellite population at higher redshift has allowed massive 
centrals to experience multiple mergers at lower redshift. 
This effect is stronger for more massive subhalos since they 
are more likely to be centrals, and they reside in higher mass 
halos with more massive satellites. 

5.2 Fraction "On" 

The fraction of halos that host recent mergers, /on , is of par- 
ticular interest for quasar or starburst evolution models, and 
quasar/starburst feedback effects such as heating of the In- 
tracluster Medium (ICM). We select mergers up to 200 Myr 
after coalescence, motivated by the expected time interval 
during which quasa rs or starbursts remain observable (e.g., 
iHopk ins et al]|2005l ). This time interval is only illustrative. 



6.1 Satellite Halo Occupation in Simulation 

The above trends are seen in Fig. [S] which shows the evolu- 
tion of the satellite occupation per halo, for satellites with 
Minf > 10"/i~^Mq in the lOO/i"^ Mpc simulation. Satellite 
occupation counts are normalized using the last output at 
z = 0.6. More massive satellites (Minf > 10^^/i"^Mq) have 
similar evolution, with a peak in the satellite halo occupation 
a.t z ~ 2.5. For a fixed satellite infall mass, less massive halos 
exhibit stronger satellite occupation evolution with redshift, 
leading to a more prominent peak. 



If we scale our /on time interval by the dynamical time, the 
slope of /on becomes slightly shallower, but the qualitative results 
do not change. 

Though if one applies a mass threshold to a population, this 
is not strictly true since mergers also scatter lower mass objects 
into the population. 
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Figure 7. Left: Fraction of halos that host a subhalo-subhalo merger within the last 200 Myr (closed points/solid lines) at z = 5 — 2.5 
for Mi„f = 10^^ — 10^'^h~^ Mq (triangles) and Mj^f = 10^^ — I0^^h~^ Mq (squares). Open points/dashed lines show the fraction of halos 
of the same mass that have suffered a major halo-halo merger in the same time interval. Right: Same, but at 2 = 1.6 — 0.6. 
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Figure 8. Evolution of the average number of satellite subhalos 
per halo for satellites with Afi^f > 10^^ Mq and several halo 
mass bins. Satellite counts are normalized to those at z = 0.6 
(the last output) and boxcar averaged across 3 outputs to re- 
duce the small number statistics noise at early times. The dotted 
black line shows fit of Eq. 119! for 3:1 mass ratio mergers using 
halo merger rate parameters from Table [T] All halos masses have 
a peak in satellite occupation at 2 ~ 2.5, and similar trends per- 
sist for higher mass satellites. Less massive halos exhibit stronger 
evolution with redshift, leading to a more prominent peak. 



6.2 Analytic Estimate of Satellite Halo 
Occupation 

The rate of change of the satellite subhalo population per 
halo, for a fixed satellite Mini, is given by the rate at which 
satellites fall into a halo (the halo merger rate) minus the 
rate at which satellites coalesce with the central subhalo 



dN, 



sat — ccn, merge 



dt dt dt 

The halo merger rate at all epochs is given by 

^-^halo, merge 



dt 



^A{l + z)°' = Aa 



(13) 



(14) 



The timescale for the satellite to coalesce with its central 



subhalo after infall is given to good approximation by 

_ Co Afhalo/Aisat 



■'sat — ccn, merge 



10 ln(l-f A/halo/Meat) 



iHubblc 



(15) 



where Co is a constant of order unity that accounts for the 
ensemble averaged satellite orbital parameters, and we leave 
it as our sole free parameter. The rate of satellite destruc- 
tion/coalescence is thus 



dJSl sat — cen, merge _ 10 ln(l -I- Mhalo/Afsat] 



dt Co Mhalo/Maat 

Combining Eqs. 1 141 and 1 161 into Eq. 1 131 one gets 



H{z) 



dN,, 



dt 



= Aa 



10 ln(l-HMhalo/Msat) 



Co 



A/halo/Ms. 

m,o 



H{z) 



(16) 



(17) 



At high redshift, H{z) 
0.1 /i Gyr~^, which implies 

diVsat _ , -a 1 ln(l-FMhalo/Msat; 



where Ho 



dt 



Aa 



Co Mhalo/Ms, 



(18) 



Across all mass and redshift regimes, the power law index for 
halo mergers is a — 2 — 2.3 (331, and we find no additional 
dependence on merger mass ratio, suggesting a universal 
power law indexlf^ Thus, the creation and destruction terms 
in Eq. [18] have differing dependencies on the scale factor, 
which implies they become equal at some redshift, where 
the satellite occupation per halo reaches a maximum. 

Using the exact evolution of the Hubble parameter (in 
our cosmology) of H{z) = Ji'o(f2m,oa~^ + ^a)^^^ in Eg.fTTl 
and integrating over a, the full evolution of the satellite oc- 
cupation per halo is 



da- 



(a + l) 



A_ 

Ho J ~" H{z 

10 Msat 



■In 



1 + ) ln(a) + K (19) 



Co Mhalo V ^^sat 

where the constant K accounts for initial/final conditions. 

Though merger rates for more discrepant mass ratios have 
higher amplitudes, as given by Eq. [5] 
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Figure [8] shows the resultant Ai'sat(a), normalized to the 
satellite occupation per halo at final output, using typical 
values for the halo merger rate from Table [l] (j4 = 0.032 
and a = 2.3). These values correspond to 3:1 mass ratio 
mergers, for satellites with Mini ~ 10^^h~^MQ this gives 
Mhaio = 3 X 1O"^"^M0. While the fit of this model to 
the simulation results is not exact, it nicely reproduces the 
general trends, especially given the simplicity of the model, 
which ignores halo/central mass growth, satellite-satellite 
mergers, and switches. For instance, it correctly produces 
a lower peak in satellite occupation for halos with satel- 
lites of more discrepant mass ratios (more massive halos for 
a fixed satellite mass, or less massive satellites for a fixed 
halo mass), relative to the amplitude at low redshift. This is 
because the decreased infall times for smaller satellite-halo 
mass ratios cause a more dramatic fall (after the peak) in 
the satellite population for lower mass halos. 

Agreement of this model with our simulations also re- 
quires Co ~ 2 ^ whi ch ca n be compared wi th other work. 
IZentner et"al] 120051 ) and Ijiang et al] (|2008l 'l find that the 
ensemble averaged satellite orbital circularity distribution 
is given by (e) — 0.5 ± 0.2, with no strong dependence on 
redshift or satellite-halo mass ratio. When applied to de- 
tailed dynamical friction timescale fits from s imulation, this 
yields Co » 0.6 jBovlan-Kolchin et al.ll2008h and Co ^ 1.4 
(|Jiang et al. |[2ooi)Er Taken at face value, our even higher 
value of Co means that our satellites are taking longer to 
merge, suggesting that the discrepancy does not arise from 
artificial over-merging in our simulations. However, exact 
comparisons are difficult given the simple nature of our an- 
alytic model, and because both of thes e groups use differ - 
ent halo and subhalo finding algorithms. I Jiang et al.l (|2008h . 
who used a simulation of roughly similar volume and mass 
resolution to ours, also incorporated hydrodynamics, which 
is likely to shorten the merger timescale since it introduces 
further dissipational effect s to the subhalo orbits and re duces 
mass loss. Compared with lBovlan-Kolchin et al] (|2008l ), who 
performed much higher resolution simulations of isolated 
halo mergers, it is possible that our satellite subhalos ex- 
perience more severe mass stripping upon infall, decreas- 
ing their mass and thus extending their subsequent infall 
time (see Eq. |5J. A more detailed investigation of satellite 
infall timescales in a cosmological setting is needed, stud- 
ies now in progress are targeting in particular the role of 
hydrodynamic eff ects (|Dolag et al.i i2008i : iSaro et al.i i2008i : 
ISimha et al.1l2008l ). 



6.3 Comparison to Other Work on Satellite 
Occupation Evolution 

Figure [S] shows a peak in the number of satellites per halo 
at z ~ 2.5. Fundamentally, the reason for this peak is that 
we select subhalos of fixed minimum infall mass in halos of 
fixed mass across time. If instead we examine the satellite 
occupation for halos above a minimum mass cut, the growth 
of massive halos (hosting more satellites) at late time would 
overwhelm the drop in the satellite population at a fixed halo 



mass, so the satellite halo occupation would grow monoton- 
ically and appear much like the satellite fraction in Fig. |3l 
which ignores halo massQ 

These results agree with the interpretation that more 
massive halos have later formation times and longer satel- 
lite inf all times, and thus host more substructure at a given 
epoch l|van den Bosch et al.ir2005l : IZentner et al.ll2005i ') . Sim- 
ilarly, for a fixed satellite infall mass, the satellite halo 
oc cupation evolves mo re ra pidly for less massiv e halos, 
as IZentner et all (120051 ) and lPiemand et all (|2007l ) found, 
though their results were based upon subhalo instantaneous 
mass and maximum circular velocity (we compare these to 
Afinf in the Appendices). 

However, the peak in satellite halo occupation in Fig. [8] 
does not appear in Halo Occupa t ion D istribution (HOD) 
evolution studies bv lZentner et al.| (l2005 l) because they use a 
fixed cut on instantaneous maximum circular velocity across 
time. As shown in Appendix A, fixed K.max probes lower 
mass at higher redshift, and this evolution in satellite mass 
overwhelms the satellite evolution of Fig. [S] leading to a 
monotonic rise in the satellit e halo occupation with red- 
shift. Similarly, IConrov et all (|2006l ) noted that the HOD 
shoulder (the halo mass where a halo hosts only a central 
galaxy) becomes shorter at higher redshift as an increas- 
ing fraction of low-mass halos host more than one galaxy, 
finding a monotonically increasing satellite population with 
redshift. However, they compare fixed satellite number den- 
sity (not mass) across redshift, which corresponds to a lower 
subhalo mass at higher redshift. Again, this overwhelms the 
evolution of Fig. [S] leading to a monotonic rise in the satel- 
lite population per halo with redshift. 

In their semi-analytic model matched to simulation, 
Ivan den Bosch et al. I l|2005h find that the average subhalo 
mass fraction of a halo always decreases with time, and they 
claim that, as a result, the timescale for subhalo mass loss 
(approximately the dynamical/infall time) is always smaller 
than the timescale of halo mass accretion (mergers). This 
implies that the satellite infall rate is always higher than the 
halo merger rate, and so the satellite HOD always decreases 
with time. However, this result refers to the total subhalo 
instantaneous mass per halo, not galaxy counts based on 
infall mass. If instead we examine the evolution of the satel- 
lite occupation per halo as in Fig. [8] selecting the satellites 
based on instantaneous subhalo mass instead of infall mass, 
we find a monotonic increase in the satellite occupation with 
no peak, in agreement with other authors above. 

These examples all illustrate how the evolution of the 
HOD is dependent both on satellite mass assignment and se- 
lection of fixed mass vs. circular velocity vs. number density 
across time. 



7 SUMMARY AND DISCUSSION 

Using high-resolution dark matter simulations in cosmolog- 
ical volumes, we have measured the rates, counts, and types 
of subhalo (galaxy) major mergers at redshift z = 0.6 — 5, 
describing their populations in terms of centrals/satellites 



■^^ The Co values of these two fits agree onl y in extre me cases, ^■^ Though the satellite halo occupation would have a higher am- 

i.e. maximally circular orbits foriBoylan-Kolchin et all 1200^) or plitude since it measures fisat/ficcntral while the satellite fraction 



maximally eccentric orbits for 0^ et al.l g measures risat/ris ubhalo — "^sat /("central +nsat). 
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and contrasting their merger properties with those of ha- 
los of the same (infall) mass. We assign subhalos their mass 
at infall (with the capacity for mass growth during satellite 
mergers) , motivated by an expected correlation with galaxy 
stellar mass, but include no further semi-analytic galaxy 
modelling. We select mergers requiring 3:1 or closer infall 
mass ratios, motivated by the expectation that these can 
trigger activity such as quasars, starbursts, and related ob- 
jects such as Lyman break galaxies, submillimeter galaxies, 
and ULIRGs. We highlight our main results as follows: 

• The merger rate per object of subhalos is always lower 
than that of halos of the same (infall) mass. Galaxies exhibit 
stronger mass dependence on the amplitude of their merger 
rate than halos, with more massive galaxies undergoing more 
mergers. While the slope of the halo merger rate per object 
is essentially redshift independent, the slope of the galaxy 
merger rate is much shallower than that of halos at z > 2.5 
and parallels that of halos at z < 1.6. 

• These differences in halos and subhalo merger rates 
arise because (1) halo mergers add mass to halos instantly, 
while central subhalo mass grows more gradually after a 
halo merger; (2) halo major mergers do not necessarily lead 
to subhalo major mergers, since central subhalos experience 
mass growth while a satellite subhalo's infall mass typically 
remains constant as it orbits; and (3) the satellite subhalo 
fraction grows with time, and satellites are twice as likely to 
be recent mergers as centrals of the same infall mass. 

• 15% — 35% of all recently merged subhalos are satellites, 
though a significant fraction of these arise from satellite- 
central parents during switches. 5% — 15% of galaxy mergers 
arise from satellite-satellite parents, with a higher fraction 
at lower redshift. 

• At 2 = 2.6 (z = 1), 30% (8%) of galaxies have experi- 
ence at least one major merger in the last 1 Gyr, regardless 
of mass. Halos are more likely to have experience multiple 
mergers in their recent history. 

• The likelihood of a halo to host a recently merged 
galaxy, /on, does not evolve with time at z > 2.5 and falls 
with time at z < 1.6. 

• Comparing galaxy and halo merger rates allows one to 
understand the evolution of the satellite halo occupation, 
and we approximated this behaviour analytically including 
fits to our simulations. Selecting subhalos on fixed infall 
mass, the satellite halo occupation for halos of a fixed mass 
increases with time at high redshift, peaks at z ~ 2.5, and 
falls with time after that. This implies similar evolution for 
the satellite galaxy component of the Halo Occupation Dis- 
tribution. 

Our results, based entirely on the dynamics of dark mat- 
ter, represent an important but preliminary step towards 
quantifying the nature of galaxy mergers in hierarchical 
structure formation. To compare to many observables we 
would need to include baryonic effects, and indeed such ef- 
fects can provide corrections to the merger rates themselves 
l|Dolag et al.ll2008l : IJiang et al.ll2008l : ISaro et al.ll2008l ). The 
merger rates here also do not include whether the subhalos 
are gas-rich (required for some observables) or not, though 
at the high redshifts we examine, we expect almost all galax- 
ies to be gas-rich. While satellites can be strippe d of much of 
their gas before mergin g with their central (e.g. JPolag et all 
Saro et aU 120081 ). we have considered primarily mas- 



sive satellites (relative to their host halos) , which have short 
infall times and thus experience less gas stripping. 

Timescales between observables and our measured 
merger event also play a role. In simulations, a quasar can 
appear up to ~ 1 Gyr after galaxy coal escence, though star- 
bursts may occur more quickly (e.g ., iHopkins et al.l l2005l : 
IStaringel et al. I l2005l : ICox et al.ll20oi ). However, morpholog- 
ical distur bance is c learest during first passage and final co- 
alescence l|Lotz et al. 2008, '). The time scales for each signa- 
ture to commence and/or persist also have a large scatter, 
ranging from 0.2 to 1.2 Gyr after the merger. Finally, spe- 
cific observations will also have specific selection functions. 
The quantitative measurements provided here provide start- 
ing points for these analyses in addition to helping to under- 
stand the properties of galaxies and their mergers in general. 

Unfortunately, our predictions are not easily compared 
to observations which estimate merger rates at 2 < 1 using 
close galaxy pair s or disturbed morpholog i es (most recently, 
Bell et al.l l2006l: iKampczvk et al] l2007l: fKartahepe et all 



20071: iLin et al.H2008l: iLotz et al.ll2008l : iMcIntosh et allboo"! 



their gi 
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Patton &: Atfield |2008| ) since they have found that the close 
pair fraction evolves as (1 -I- z)" with a diverse range of ex- 
ponents from a = to 4. Furthermore, translating these 
observations into galaxy or halo merger rates requires in- 
cluding more physical effects, e.g. time dependent galaxy 
coalescence timescales l|Kitzbichler fc White! l2008l : iMateud 
12008' ) or inclusion of changes in numbers of host halos with 
redshift l|Berrier et al.ll2006l) . 

While we were p reparing this work for prepara- 
tion, ISimha et al.l (|2008l ) appeared which considers detailed 
merger properti es of satellite s ubha los in a hydrodynamic 
simulation, and lAngulo et al.l l|2008l ) appeared which con- 
siders satellite mergers in t he Millennium simulation. After 
this paper was submitted, IStewart et al.l l|2008f ) appeared, 
which also considers merger rates as a function of mass, 
redshift, and mass ratio. 
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APPENDIX A: SUBHALO MASS AND 
CIRCULAR VELOCITY 

Here we further elaborate on the details of subhalos in our 
simulations, including their radial density and circular ve- 
locity profiles and the relation between mass and maximal 
circular velocity, including its evolution with time. 

Figure lAll shows the radial density and circular ve- 
locity profiles for 8 satellite and 8 central subhalos with 
K.max =i 250 km/s (Af ~ lO^^fe-^Afo) at 2 = 2.6 in 
the 100/i~^ Mpc simulation. Circular velocity is defined as 
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Figure Al. Left: Radial density profiles (top) and circular velocity profiles (bottom) of all matter around 8 satellite (dotted) and 8 
central (dashed) subhalos with 14, max — 250 km/s at ^ = 2.6. Right: Same, but only for matter assigned to the subhalos. 
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Figure A2. Relation between subhalo infall maximum circular 
velocity, infi and subhalo infall mass, Mi^f, at z = 2.6. Black 
points show a 25% sub-sample of all subhalos as a measure of 
scatter, the solid red line shows the least squares fit to Eq. lAll 
and the dashed red lines show the Icr scatter, for subhalos with 
Afinf > lO"/i-iM0. 



Vc = y/GM{< r)/r, and since subhalos follow NFW density 
profiles, with a break in the power-law density profile at the 
scale radius, r^, they have a maximum value in their circu- 
lar velocity profiles, 14, max, at rmax = 2.2rs. The left panels 
show the radial profiles of all matter surrounding the subha- 
los, while the right panels show only that of matter assigned 
to the subhalos. Satellites and centrals have similar profiles 
at small radii, though the right top panel shows that satel- 
lites exhibit signs of tidal truncation at ~ 50/i~^ kpc. This 
is also visible in the circular velocity profile in the bottom 
left panel, where Vc for satellites rises sharply, exhibiting a 
transition to their host halos. For the centrals, the rise in Vc 
beyond ~ lh~^ Mpc arises from neighbouring structures. 
Figure IA2I shows the relation between subhalo infall 
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Figure A3. Evolution of the ratio of subhalo infall maximum 
circular velocity to infall mass, B{z) = V^^^f/Mi^l (solid), and Icr 
scatter (dashed) for subhalos with Mj^f > 10^^/i~^Mq. Dotted 
line shows fit to B{z) of Ea.lA2] 

mass, Minf, and infall maximum circular velocity, Vdni, for 
subhalos at z = 2.6. We fit this relation to 

KTmf = B(^)Mnf (Al) 

for all subhalos above 10^^ Mq, finding 7 = 3 holds to 
good approximation at all redshifts we examine, in agree- 
ment with the virial relation Vi?niax oc M/R oc M/M^^^ oc 
M^^^. The outliers with large M-mi relative to T4,inf are satel- 
lites that experienced a major merger; under our prescrip- 
tion, a satellite child's M^^f is the sum of its parents' Minf, 
but a child's Vc,inf is that of its highest Vc,ini parent. 

Fixing 7 = 3, Fig. IA3I shows the evolution of the am- 
plitude B{z) = V^^inf/Minf. A subhalo of a given mass has 
a higher maximum circular velocity at higher redshift, re- 
fiective of the increased density of the universe when the 
subhalo formed. As a subhalo subsequently accretes mass. 
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Figure Bl. Satellite subhalo (scaled) mass function, for various 
host halo mass bins, at ^ = 1. Thick curves show satellite mass 
selected on Mj^f, while thin curves show satellite mass selected 
on instantaneous bound mass. While the instantaneous bound 
mass function exhibits no dependence on halo mass, the infall 
mass function has a higher amplitude for more massive halos. 
An appreciable number of satellites exists at Msat,inf/A^halo ~ 1 
because of switches. Below the rollover at high satellite mass, both 



mass functions scale as 



,, , , f — =7t-r r (X M^Y'" in agreement 

with instantaneous satellite subhalo bound mass functions found 
by numerous authors (see text). 



its 14 max grows Hiore slowly than it mass (in cases of slow 
mass growth, we fi nd Vc.msoc can r e main constant). This is 
in agreement with IPiemand at al.l (|2007h , who found that 
halos undergoing mild mass growth (no major mergers) had 
less than 10% change in K;,max and rmax- We find that the 
evolution of B{z) (within the redshifts we probe) can be 
well-approximated by 

B{z) = 6.56 X 10~''e°-^''"(km/s)^/iMQ\ (A2) 

Thus, Minf = 10" (lO^^)/i"^M0 subhalos correspond 
to K.inf =i 120 (250) km/s a.t z = 2.6 and K,inf 
100 (200) km/s a.t z — 1. Conversely, for fixed V^.inf, a sub- 
halo is about half as massive ai z — 2.6 than at z — 1. 

Since we fit relations of satellite subhalo properties at 
infall, our results above are applicable equally to satellite 
and central subhalos. In addition, our results change by only 
a few percent if we instead consider host halos. 



APPENDIX B: 
FUNCTION 



SATELLITE SUBHALO MASS 



Figure iBll shows the satellite subhalo (scaled) mass function 
vs. the ratio of satellite mass to halo mass, for various host 
halo mass bins ai z = 1. Thick curves show satellite masses 
selected on Minf, while thin curves show satellite masses 
selected on instantaneous bound mass. The rollover in the 
Mini curves at low satellite mass indicates where satellites 
become numerically disrupted by resolution effects. Numer- 
ical disruption occurs at higher satellite Minf for more mas- 
sive halos, indicating that satellites of a fixed Minf experi- 
ence more pronounced tidal stripping in higher mass halos, 
where dynamical friction timescales and central densities are 
higher. The rollover in the highest halo mass bin occurs at 



Msat.inf ~ lO^^/i~^M0, which sets our minimum subhalo 
mass for robust tracking. 

We find that the scaled instantaneous bound mass func- 
tion exhibits little-to-no systematic depend ence on halo 
mass, in agreement with lAngulo et all (|2008hF^ This is in 
contrast to the scaled infall mass function, which shows more 
satellites at a given mass ratio for more massive halos. This 
difference is driven by subhalo mass stripping. Averaged over 
the entire satellite population at this redshift, the satellite 
instantaneous bound mass is ~ 30% that of Minf, as can 
be seen by the x-axis offset of the solid and dashed curves. 
However, satellites exhibit less average mass loss in low-mass 
halos than high-mass halos, the instantaneous to infall mass 
ratios being 40% and 25%, respectively. This is consider- 
ably higher than the 5% — 10% at z = 1 foun d in the semi- 
analytic model of Ivan den Bosch e t al.' ('2005') . Additionally 
van den Bosch et al. (20051) andlCiocoh ct al. (2008) found 
the opposite trend with halo mass, i.e., that the average 
mass loss of satellites is higher for lower mass halos. They 
find that this arises because lower mass halos form (and ac- 
crete their subhalos) earlier, when the dynamical/stripping 
timescale is shorter. Thus, the satellites of lower mass halo 
are strip ped both more rapidl y and over a longer time period 
(see also lZentner et ahllioosl ). Finally, ICiocoli et J] (|2008| ) 
find that the scaled mass functions of subhalos at infall does 
not depend on halo mass, in seeming contrast with Fig. IBll 

These discrepancies likely arises because we examine 
the masses of extant sub halos in our simu l ation, while 
Ivan den Bosch et all (|2005l ) and iGiocoh et al] (|2008l ) track 
subhalo mass loss much longer than our simulation does. 
Their semi-analytical model of subhalo mass loss has no pre- 
scription for central-satellite mergers, which preferentially 
serve to reduce highly stripped satellites from our sampler'*! 
Since satellites of a given infall mass to halo mass ratio have 
lower mass in lower mass halos, they are closer to our mini- 
mum subhalo finding mass threshold. Thus, a fixed amount 
of stripping will cause lower mass halos to have a reduced 
population of satellites of a given infall mass to halo mass 
ratio. It is unclear at what level of subhalo mass stripping 
we should expect the galaxies they host to become disrupted 
as well. 

Various studies using high resolution simulations have 
explored in detail the slope of the subhalo mass function 
("De Lucia et al."2004'; 'C ao et al.ll2qo3: iDiemand et alllioOTl : 
.Madau et al.. ^0&; .Angulo et all |2008| ). finding that the 
(instantaneous) mass function of subhalos goes as A'sat(> 
Meat) oc (Mhaio/M,at)", with a = 0.9 - iQ We note, 
though, that these fits are based on cuts on instantaneous 
subhalo mass or circular velocity, not those at infall. How- 
ever, we find that, between the rollover in Afinf at low and 
high satellite mass, the slopes of the mass functions selected 
on infall and instantaneous bound mass are the same within 
error, and a — 0.9 fits our mass function selected on ei- 
ther mass. The similarity of the two slopes also implies that 



We thank the refere e for su gges ting we show this qu antity. 
24 See iTavlor fc Babul] ||2005|) and IZentner et all (|2005|) for de- 
tailed comparisons of simulated subhalos with analytic models. 
2^ Afaat(> Msat) has the same power law dependence on subhalo 

"^^^ dln(Maat/Mhalo)) ■ 
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the mass loss rate of satellites do es not depend stron gly on 
satellite mass, in agreement with lGiocoli et al] (|2008l ). 
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